Active transport in dense diffusive single-file systems 
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We study a minimal model of active transport in crowded single-file environments which gener- 
alises the emblematic model of single file diffusion to the case when the tracer particle (TP) performs 
either an autonomous directed motion or is biased by an external force, while all other particles 
of the environment (bath) perform unbiased diffusions. We derive explicit expressions, valid in the 
limit of high density of bath particles, of the full distribution P n (X) of the TP position and of all 
its cumulants, for arbitrary values of the bias / and for any time n. Our analysis reveals striking 
features, such as the anomalous scaling oc y/ri of all cumulants, the equality of cumulants of the 
same parity characteristic of a Skellam distribution and a convergence to a Gaussian distribution 
in spite of asymmetric density profiles of bath particles. Altogether, our results provide the full 
statistics of the TP position, and set the basis for a refined analysis of real trajectories of active 
particles in crowded single-file environments. 
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Introduction. Single-file diffusion refers to one- 
dimensional diffusion of interacting particles that can not 
by-pass each other. Clearly, in such a geometry the ini- 
tial order of particles remains the same over time, and 
this very circumstance appears so crucial that the move- 
ments of individual particles become strongly correlated: 
the displacement of any given tracer particle (TP) on 
progressively larger distances necessitates the motion of 
more and more other particles in the same direction. 
This results in a subdiffusive growth of the TP mean- 
square displacement X 2 ~ yi, first discovered analyt- 
ically by Harris [T] and subsequently re-established for 
systems with differently organised dynamics (see, e.g., 
Refs.|2 -8 ). Nowadays, a single-file diffusion, prevalent 
in many physical, chemical and biological processes and 
experimentally evidenced by passive microrheology in ze- 
olites, transport of confined colloidal particles or charged 
spheres in circular channels [9 13 . It provides a paradig- 
matic example of anomalous diffusion in crowded equilib- 
rium systems, which emerges due to a cooperative many- 
particle behavior. 

On the other hand, systems that consume energy for 
propulsion - active particle systems - have received grow- 
ing attention in the last decade, both because of the new 
physical phenomena that they display and their wide 
range of applications. Examples include self-propelled 
particles such as molecular motors or motile living cells 
[2] , and externally driven particles, such as probes in ac- 
tive microrheology experiments |15| . The intrinsic out-of- 
equilibrium nature of these systems leads to remarkable 
effects such as non-Boltzmann distributions [TB], long- 
range order even in low spatial dimensions 17 and spon- 
taneous flows [18]. In particular, ID assemblies of active 
particles have been extensively studied in the context of 
TASEP models. 



However, up to now, active transport in diffusive sin- 
gle file systems, which involve an active TP performing 
an autonomous directed motion or pulled by a constant 
external force / in a ID bath of unbiased diffusive par- 
ticles with hard-core interactions, has drawn uncompa- 
rably less attention |29j . Such dynamics, depicted in 
figure [T] provides a minimal model of active transport 
in crowded single-file environments, which schematically 
mimics situations as varied as the active transport of a 
vesicle in a crowded axone [19], directed cellular move- 
ments in crowded channels |20j or active-microrheology 
in capillaries [T5]. In this context, the only available the- 
oretical results concern the large time behavior of the 
mean displacement X of the TP, which has been shown to 
grow sub-linearly with time X ~ \/t j2"D - f!M] . In fact, the 
biased TP drives the bath particles to a non- equilibrium 
state with an asymmetric distribution: the bath parti- 
cles accumulate in front of the TP thus increasing the 
frictional force, and are depleted behind. The extent of 
these perturbations grows in time in proportion to \fi 
and characterizes a subtle interplay between the bias, for- 
mation of non-equilibrium density profiles and backflow 
effects of the medium on the TP. 

In this Letter we focus on this minimal model of active 
transport in diffusive single file systems. Going beyond 
the previous analysis of the TP mean displacement, we 
present in the limit of high density of bath particles ex- 
act expressions of the full distribution P n (X) of the TP 
position and of all its cumulants, for arbitrary values of 
the bias / and for any time n. In particular, in this 
high density limit, it is shown that at large times this 
distribution is Gaussian, with mean X ~ af{p)^/n and 
variance growing asymptotically as *Jn. Remarkably, in 
this limit the variance is proved to be independent of /. 
Altogether, our results provide the full statistics of the 
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TP position, and set the basis for a refined analysis of 
real trajectories of active particles in crowded single-file 
environments. 

The model. Consider a one-dimensional, infinite in 
both directions line of integers x, populated by hard-core 
particles present at mean density p, performing symmet- 
ric random walks. At t = we introduce at the origin 
of the lattice an active TP, hopping on its right (resp. 
left) neighbor site with probability p\ (resp. P-i), which 
process is also constrained by hard-core exclusion. In 
what follows, we focus on the limit of a dense system, 
corresponding to the limit of a small vacancies density 
po = 1 — p <C 1. In this limit it is most convenient 
to follow the vacancies, rather than the particles. We 
thus formulate directly the dynamics of the vacancies, 
which unambiguously defines the full dynamics of the 
system. Following 25, 26J, we assume that at each time 
step each vacancy is moved to one of its nearest neigh- 
bours sites, with equal probability. As long as a vacancy 
is surrounded only by bath particles, it thus performs a 
symmetrical nearest neighbors random walk. However, 
due to the biased nature of the movement of the TP, 
specific rules have to be defined when a vacancy is adja- 
cent to the TP. In this case, if the vacancy occupies the 
site to the right (resp. to the left) of the TP, we stip- 
ulate that it has a probability q\ = l/(2p 1 + 1) (resp. 
g_i = l/(2p_i + 1)) to jump to the right (resp. to the 
left) and 1 — qi (resp. 1 — q-i) to jump to the left (resp. 
to the right). These rules are the discrete counterpart of 
a continuous time version of the model [3D], as shown in 
[25] , Note that a complete description of the dynamics 
would requires additional rules for cases where two vacan- 
cies are adjacent or have common neighbours; however, 
these cases contribute only to 0(p 2 ), and can thus be left 
unstated. 



Following Refs. 



], we then represent pV^(X) as: 
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FIG. 1: (color online) Model notations. 



Single file with a single vacancy. We start with an aux- 
iliary problem in which the system contains just a single 
vacancy initially at position Z, and which will be proved 
next to be a key step in the resolution of the complete 
problem with a (small) concentration of vacancies. Let 
p r z(X) denote the probability of having the TP at site X 
at time moment n, given that the vacancy commenced 
its random walk at Z. Clearly, in a single vacancy case 
this probability is not equal to zero only for X = and 
X = 1, if Z > 0, and X = and X = -1, if Z < 0. 



p n z (X) = S x , 

V i=o 

p — 1 7TOl,7Ti.2s" . ,mp=l m p -j-i=0 
rn p+1 \ 

1 E F (-iy J xF™» 1)P+1 ...F™F™\ (1) 

where 5 a j, = 1 when a = b and is equal to zero, otherwise, 
and F£ is the probability that the vacancy, which started 
its random walk at site Z, arrived to the origin for the 
first time at time moment n. 

Let now <?(£) denote the generating function of any 
time-dependent function g(n), <?(£) = Y^Lodi 71 )^"- 
Then, Eq. ([I]) implies that the generating function of 
the propagator of the single-vacancy model can be ex- 
pressed via the generating functions of the corresponding 
first-passage distributions as 

, y n S x , (l-F ±1 ) + 6x,±iF ±1 (l-F Tl ) 

where we have used the short notations F±i = F±i(£). 

Single file with a small concentration of vacancies. We 
now turn to the original problem with a small but finite 
density po of vacancies and aim to express the desired 
probability P n (X) of finding the TP at site X at time 
n via the propagator for a single-vacancy problem. We 
consider first a finite chain with L sites, M of which are 
vacant and the initial positions of the latter are denoted 
by Zj, j = 1, . . . , M. Then, the probability P n (X\{Zj}) 
of finding the TP at position X at time moment n as a 
result of its interaction with all the vacancies collectively, 
for their fixed initial configuration, writes 

PniXliZ,}) = Yl 5 X y l + ... +YM P n ({Y J }\{Z ] }), 
Y U Y2,...,Y M 

(3) 

where P n ({Yj}\{Zj}) stands for the conditional probabil- 
ity that within the time interval n the TP has performed 
a displacement Yj due to interactions with the first va- 
cancy, a displacement Y% due to the interactions with the 
second vacancy, and etc. In the lowest order in the den- 
sity of vacancies, the vacancies contribute independently 
to the total displacement of the tracer, so that the latter 
conditional probability decomposes 



M 

Po->0 * A J 



(4) 



where O^j ) ^ s the single- vacancy propagator and the 
symbol ~ signifies the leading behavior in the small 

Po-s-0 

density of vacancies limit. Note that such an approxi- 
mation yields results which are exact to the order O(po), 
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and hence, such a description is expected to be quite ac- 
curate when po <C 1 [25, 26J. Next, we suppose that 
initially the vacancies are uniformly distributed on the 
chain (except for the origin, which is occupied by the TP) 
and average P n (X\{Zj}) over the initial distribution of 
the vacancies. In doing so and subsequently turning to 
the thermodynamic limit, i.e. setting L oo, M — > oo 
with M/L = po kept fixed, we find that the generating 
function of the second characteristic function 



(5) 



active transport in dense diffusive single file systems as 
detailed below. 

(i) First conclusion we can draw from Eq. (11) is that 
for arbitrary / (including / — 0) all odd cumulants have 
the same generating function ft dd(£)i an d a H even cumu- 
lants have the same generating function /c C vcn(£)- This 
means that at any moment of time and for any / all cu- 
mulants Ktj\ (n) with arbitrary odd j are equal to each 
other, K(2j+\){n) = K dd(^), and so do all the cumulants 
with arbitrary even j, K( 2 j)i n ) = K CVO n(^). 

Parenthetically, we note that, in the classical case of 
single file diffusion (i.e. / = 0), the generating function 
in Eq.(ll ) can be inverted explicitly to give n odd (n) = 



satisfies 
lim^Ml 



and for arbitrary time moment n 



E 

£=±1 



z=i 



(6) 

Our last step consists in the explicit determination of F±i 
and Ez=i m Eq.M. We note that both can be 

readily expressed via the nrst-passage time density at the 
origin at time n of a symmetric one dimensional Polya 
random walk, starting at time at position I, denoted 
as since, by partitioning over the first time when the 
sites adjacent to the origin are reached, we have: 



q±i)S n .i +q±i^2f k ~ 

k=l 



\1)F__ 



±i 



(7) 



Multiplying both sides of Eq. Q by £ n , performing 
summation over n and taking into account that /;(£) = 
E /PC = (C 1 - V 1 -C 2 )/?) 1 ' 1 we find that 



F±i = 



(1 - q±i)Z 



Similarly, noticing that 



b z 



ELi/I-i^r fe if^>o, 
ELi^i-^"r fc if^<o, 



(8) 



(9) 



and using the definition of /;(£) given above, we obtain 



(10) 



Gathering the results in Eqs. ([7| to ( |10| ), substituting 
them into Eq. ([6| , we finally derive our central analytical 
result which defines the exact (in the leading in po order) 
generating function of the cumulants of arbitrary order 



«(j)(0 



F 1 (l-F_ 1 ) + (-lp'F_ 1 (l-F 1 ) 



lim 

Po^O po (l_£)(l_(l_0-Z£2)/£)(l_ Fl J7_ 1 ) 

(11) 

This result gives access to the full statistics of the position 
of the TP and puts forward striking characteristics of 



lim ■ 

Po^O 



2 r(L^i] + f) 
- V5Fr(|*^J + i) 



(12) 



where T(-) is the Gamma function and [a;J is the floor 
function. This expression, which can be shown to be 
compatible with the well-known Gaussian form in the 
large time limit, seems to be new. 

(ii) Second, turning to the limit £ — > 1 (large-n limit) 
we find the leading in time asymptotic behavior of the 
cumulants of arbitrary order: 



(2j+l) 



lim 

Po^O po 



lim 



Oi -P-x)\j ~ -2piP-i(Pi ~P-i) + o(l) 



Po->0 po 



+ o(l), j = 0,1,2, 



(13) 



(14) 



Equations (13 1 and (14) signify that, remarkably, the 



leading in time behavior of all even cumulants is inde- 
pendent of the force /, while the leading in time behavior 
of all odd cumulants does depend on /. In addition, for 
the standard choice of the transition probabilities such 
that pi = 1 — P-i and p\jp-\ = exp(/3/), where j3 is the 
reciprocal temperature, and for the specific case j = 0, 
we check from Eq. ( 13 ) that 



X 



lim — = tanh(/3//2) x /2n/7r, 

Po^O po 



(15) 



which reproduces, for j — 0, the results of [22] and [23) . 
Note that this anomalous scaling oc y/n holds for all cu- 
mulants. 

(iii) We finally provide an explicit expression of the 
full distribution function P n (X) for any n. As a matter 
of fact, the equality at leading order in po of cumulants 
of the same parity proved in point (i) shows that the 
distribution associated to these cumulants is of Skellam 
type [28!, so that : 



Pn{X) ~ 

fir, 



~ exDf- K (n)) ( *even(n) + K odd (n) \ X/2 

— ex P \ K cvcn [n) ) i - - - r I 

Po^O V K ovon(^) ^ K odd (n)J 

X Ix[sl «2ven(»») - K odd( n ) J ' ( i6 ) 
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where Ix{~) is the modified Bessel function. Importantly, 
we find that despite the known asymmetry of the con- 
centration profile of the bath particles [53] , the rescaled 
variable (X n — K dd{n))/^/ K cvcn (n) is asymptotically dis- 
tributed accordingly to a normal law. More precisely, the 
convergence to this Gaussian distribution can be quan- 
tified by the skewness 71 = K(3){n)/K%\{ri) and excess 
kurtosis 72 = K^(n)/ K^(n) of the distribution P n (X). 
From Eqs.(13) and (14) we readily find that in the lead- 



ing in po order 

tanh(/S//2) 



(JL)"* + „,!/„./.,, (17) 



71 



1/2 



and 



72 



Po- 



1 

>0 p 



1/2 



+ o{l/n 1 ' 2 ). 



(18) 



Note that 71 > which signifies that the right tail 
of P n (X) is longer and the fluctuations are more pro- 
nounced for X > X where the bath particles accumu- 
late, than in the region X < X depleted with the bath 
particles. 
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FIG. 2: (color online) Odd cumulants at time n = 100 vs 
po • The straight lines define our predictions in Eq. (Ill for 
different values of pi , 



while the filled and empty symbols are 
the results of numerical simulations for the first and third 
cumulants, respectively. Circles are results for p\ — 0.55, 
squares - for pi = 0.6, triangles - for pi — 0.75 and diamonds 
- for pi = 0.98. The inset shows analogous results for the 
second and the fourth cumulants. 

Note finally that the regime of validity of our expres- 
sions with respect to the density po is tested in Fig. [2] 
where we compare our theoretical predictions for the cu- 
mulants, obtained by the inversion of our general Eq. 



( 11 ), against the results of numerical simulations for dif- 
ferent values of the density po of the vacancies, for dif- 
ferent forces / (defined as j3f — ln(p 1 /p_ 1 ) and a fixed 
time moment n = 100. We observe a very good agree- 
ment for very small values of po and conclude that, in 
general, the approach developed here provides a very ac- 
curate description of the TP dynamics for p < 0.1. Fur- 
ther on, in Fig. [3] we plot our theoretical predictions for 



the time-evolution of the cumulants for different values of 
the force and at a fixed density po. Again, we observe a 
perfect agreement between theory and simulations. Note 
that for small fields the reduced odd cumulants approach 
1 from above, while for strong fields from below. Last, we 
compare in Figj4] our prediction in Eq. ( 16 1 against the 
numerical data and again observe a very good agreement 
between our analytical result and numerical simulations. 
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FIG. 3: (color online) Reduced cumulants K( cvcn )(n) = 
W( CV cn) (w)/-y /2n/7r and K (odd) (n) = K (od d) (™)/[(pi - 
P-i)y / 2n/ir — 2pip-i(pi — p-i)] vs time n for po = 0.01 and 
a) pi = 0.6 and b) pi = 0.98. Solid lines give the results of 
the inversion of Eq. ( |11[ ), while symbols are the results of 
numerical simulations. 



• n = 10 4 

• n = 10 5 

a n = 10 6 
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FIG. 4: The distribution P n (X) for p = 0.002 and pi = 0.98 
The dashed lines are our theoretical predictions in Eq. ( 16 1 
while the symbols are the results of numerical simulations. 



Summary. We have solved a minimal model of ac- 
tive transport in crowded single-file environments. Our 
approach generalises the emblematic model of single file 
diffusion to the case of an active TP. We have derived 
explicit expressions, valid in the limit of high density of 
bath particles, of the full distribution P n (X) of the TP 
position and of all its cumulants, for arbitrary values of 
the bias / and for any time n. Our analysis reveals strik- 
ing features, such as the anomalous scaling cx ^/n of all 
cumulants, the equality of cumulants of same parity char- 
acteristic of a Skellam distribution and a convergence to a 
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Gaussian distribution in spite of asymmetric density pro- 
files bath particles. Altogether, our results provide the 
full statistics of the TP position, and set the basis for a 
refined analysis of real trajectories of active particles in 
crowded single-file environments. 
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